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1. Introduction 

Electromagnetic (EM) properties of hadrons offer a rich source of interesting and important 
phenomena. The patterns of the mass splittings between charged and neutral mesons or the mass 
splittings among the octet or decuplet baryons are sensitive to the isospin breaking of up and down 
quark masses and the presence of EM interactions. It is also known that the width difference of 
the p + and p° mesons and the hadronic light-by-light scattering amplitude, which may not be 
computed from measured EM properties of hadrons, play an important role in the standard model 
(SM) prediction of the anomalous magnetic moment of muon. 

Recent developments in lattice QCD in both hardware and software have advanced the field 
closer to the goal of QCD calculations without approximation; large-scale, high precision un- 
quenched simulations are becoming available [JlJ]. What is important here is that the statistical 
error of pseudoscalar meson masses are well under control, with size typically a few percent or 
less. Remembering that charged-neutral meson splittings are of 0(a em ) ~ 0(1%), where a em is 
the fine structure constant, it is expected that once EM interactions are successfully included, it will 
be possible to determine the up and down quark masses, which are the poorly known parameters of 
the SM, from first principles by using such splittings as inputs, and thus we can check the simplest 
solution to the strong CP problem, m u =0. 

In this talk, we focus on the determination of the light quark masses. The strategy is basically 
following Refs. [g], where the EM fields are introduced in a non-compact form, and combined 
with QCD link variables to realize QCD + QED theory on the lattice. While this pioneering work 
made use of the Wilson quark action in the quenched approximation to QCD, here we employ 
domain wall fermions [|3|, and QCD configurations with two flavors of dynamical domain wall 
fermions, generated by the RBC Collaboration Our activity toward the lattice calculation of 
the anomalous magnetic moment of muon is reported in Ref. [0]. 

2. QCD+QED system 

We employ unquenched QCD gauge configurations t/ qc d,/x(*)> generated with parameters V = 
16 3 x 32, L s = 12, Ms = 1.8 and three different sea quark masses, m sea = 0.02, 0.03, 0.04 [§]. 
The lattice spacing is determined to be 1.691(53) GeV from ra p =770 MeV, and hence the physical 
volume is roughly (1.9 fm) 3 . Calculations are performed on about 200 configurations from the 
5,000 trajectories available at each value of the sea quark mass. 

A non-compact form is adopted for the action of the U(l) gauge fields, A em ^(x) To gen- 
erate A em ,/i(jc), we first rewrite the action in momentum space, and impose the Coulomb gauge 
condition. After diagonalizing the kernel, the gauge fields are chosen randomly according to a 
Gaussian distribution and with e=l. The configuration A em ^(x) is then obtained by Fourier trans- 
formation. It is worth noting that with the use of the non-compact action and the generation pro- 
cedure described above there is no auto-correlation among the configurations, even for arbitrarily 
small coupling. Since quenched QED is a free theory, the fine structure constant a em receives no 
renormalization. Using A em ^(x) and the quark electric charge Q q , we obtain a t/(l) link variable, 
(U em ^(x)) Q i = e'GAm.nW. Configurations for the (QCD + QED) theory are then constructed by 
^qcd,ji(-x) x (U em n(x)) ® q , which are used in the inversion of the Dirac operator for valence quarks. 
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Though e-\ is adopted in the configuration generation, any value of the quark's electric charge can 
be realized by tuning Q q . 

In this system, the flavor non-singlet axial Ward-Takahashi (WT) identity for two flavors of 
domain wall fermions with mass m q \ and m q2 and charge Q q \ and Q q 2 is given by 



T° T 3 



= Kl + m q2) P a ^) + {m ql -m q2 ) q{x)Y 5 { — ,—}q{x)+2Jl q {x)+Y,£{s)K(x)i 

(2.1) 

where (x) takes the same form as the conserved axial- vector current in the pure QCD case, but 
with link variables, t/ qc d,/i(*)x diag.[(£/ em ,^ {x)) Qqi , (f/ em ,/j ( X )) G?2 ]» P a (x) a pseudoscalar density, 
/f ? (x) the same as in pure QCD, and 
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2 > 
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X a T 3 



, (2.2) 



and £(5) = 1 for 1 < s < L s /2 and —1 for L s /2 < s < L s . Notice that 3£ s a (x) vanishes when a = 3 
or Q q \ = Q q2 . 

From the analogy to the fact that in pure QCD with the domain-wall formalism the presence 
of the singlet axial anomaly becomes clear by considering the disconnected diagram involving 
the singlet, mid-point current, it is inferred that the similar understanding is possible for the £7(1) 
anomaly, but this time the disconnected diagram involving the flavor non-singlet, mid-point current 
J$ q (x) plays an important role. Due to the presence of the t/(l) anomaly, none of the components of 
s#n{x) is conserved for Q q \ / Q q2 even in the quark massless limit, and hence there is no Nambu- 
Goldstone (NG) boson in this system If you consider the lattice measurement of m n a through 
the standard method using two-point correlation functions, the possible origin of the non- vanishing 
rrijfi is the 7T°-singlet mixing, and it could affect m n o by 0(a^ m ) or higher. At present, we do 
not include the 7T°-singlet mixing. Thus, n° must become massless in the quark massless limit. 
Notice that even after including the mixing effects, it is justified to regard n° as an NG boson up 
to 0{al m ). In addition to the above approximation, we are neglecting the disconnected diagram 
in 7T -7r° correlation function (for now). It should be noted that the latter approximation does not 
affect the fact that "7T°" becomes massless in the quark massless limit, and more importantly, the 
point where m n a vanishes. 



3. Numerical results 



At present, we only have meson data in which the two valence quarks are degenerate with 
the sea quark mass. The electric charge of the valence quarks is set to Q q = +2e/3 for the up 
quark and — e/3 for the down quark. In this work, we take three values of e, which correspond 
to 0£e m = oc e m,phys> ^fj > . to examine the a em dependence of the meson mass splittings. The 
mass of the pions is extracted from two-point correlation functions as usual. We employ two 
different operators, P a (x) and Aq(x), as the interpolating operator. Since the EM interactions violate 
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Figure 1: The a sm dependence of the splittings. The results from (AqAq) are shown. 



isospin symmetry, 71° mixes with the lightest isosinglet state. However, this mixing turns out to start 
only at 0(a 2 m ), and hence we neglect this effect as the present statistics cannot resolve such a small 
effect. Although the calculation of m^o requires a disconnected diagram, we do not include it, as 
mentioned earlier. We will come back to this point later. 

At the leading order of chiral perturbation theory (ChPT) including the EM interactions, m 2 n+ 
and m 2 ^ are given by 

m 2 K+ = a em A+ + 2 (S + a em Ay) (m/ )+ + m res ) , (3.1) 
m 2 n a = 2 (B + a em A m) ) (m f}+ + m res ), (3.2) 



to 0(a 2 m ), where mf t ± = (m/^ ±m/ iM )/2. In the following, these functional forms are used to 
determine the low energy constants (LEC's), Bo, Ay, A+ and Aq . m res is estimated by finding 
m.f at which m 2 ^ vanishes, since n° can be considered as an NG boson. 

Before discussing the determination of the LEC's, let us mention 0(a em ) 2 corrections. Fig- 
ure [I] shows the a em dependence of the n + -n and 7i Q -n^ mass splittings, where n@ denotes the 
pion in pure QCD. It is seen that each splitting is well described by a linear function of a em in the 
range i < a em < 0.08. This means that 0(a 2 m ) contributions are not significant in this range of 
Oem- This observation allows us to treat the EM interactions perturbatively. 

Now we proceed to the determination of the LEC's. The entire analysis is done under a jack- 
knife procedure. Bo is obtained by fitting the pure-QCD pion data, m 2 ^. By fitting m 2 l+ — m^p as 
a function of mj, the combination (A^ + A^) and A^ are extracted. Finally A^ is obtained 
from the slope of m 2 n0 — m 2 ^ . In Fig. ||[ the m f dependence of m 2 n+ — rrv\ and m 2 ^ — m 2 ^ and the 
resulting fit curves are shown as examples. From these fits we obtain 

Bo = 2.09(4), A^ =0.0429(71), A? 5 = 2.28(31), A[ m) = 1.60(28), fromA , (3.3) 
#0 = 2.06(4), A^ =0.0466(51), A^ = 2.56(29), aJ' = 1.94(27), from/ 5 , (3.4) 

in lattice units. Since we have not measured masses of non-degenerate mesons, we simply apply 
eqs. ( fO| ) and (3_2) to estimate the kaon masses, 

m\+ = OemA^-l- (B + a era A^) (m /rV + m reSjS + m /jM + m reS;U ), (3.5) 
«4o = (Bo + a em Ajp ) (rrifj + m reSjS + m f j + m reS;d ) • (3-6) 
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Figure 2: The mj dependence of n + -7t° splitting (left) and 7t°-7t@ splitting (right). 



Using the lowest order ChPT expressions, 

mlo = 2 (#0 + Oem A m) ) (m /;+ + m reS!+ ) , 

(m K + 2 + m K o 2 ) = 2B Q (m /rS + m f& + m ieSiS + m reSi+ ) 



+« em + (aJ"> + A« ) {m f , + m res ,.)) . 



(3.7) 



(3.8) 



(m|o-m| + )-(m^o 



-«?„ 



2S (m/_ +m res _) + a em (aJ m) - A^) (m /:J + m res ,,),(3.9) 



the results of eqs. ( p.3| ) and ( p.4\ >, and the experimental values of the meson masses, we can deter- 
mine nif ± andnifs (0(a em nif t ±) terms have been neglected). Using a _1 =1.691 GeV, and the non- 
perturbative value of the quark mass renormalization factor 1/Z m = Z s =0.62 [§] (m MS (2GeV) = 
Z m (my +wi res )), we obtain the following preliminary results for the quark masses (statistical error 
only) 



m^(2GeV,ref) = 2.84(23) MeV (A ), 2.89(8) MeV (P), 
m^p(2GeV,ref) = 5.41(24) MeV (A ), 5.49(13) MeV (P), 
mp(2GeV,ref) = 106.8(8) MeV (A ), 106.5(5) MeV (P). 



(3.10) 
(3-11) 
(3.12) 



While the statistical errors appear to be under control, we emphasize that our results are determined 
in the quenched approximation of QED, and we have used only mj = m va \ = m sea data points in our 
fits ('ref ' in the above results is stressing it). It has been reported that in the pure QCD calculation 
partially-quenched (PQ) valence quark mass effects at the next-to-leading order in ChPT tend to 
increase the strange quark mass compared to the lowest-order extrapolation used here [§]. In ad- 
dition, the naive definition of m res used in this analysis amounts to a constant shift downward of a 
few percent on each of the quark masses compared to the results in Thus our results should be 
considered preliminary. On the other hand, the effects of the EM interactions relative to pure QCD 
(and on the m u — mass difference itself) should not depend significantly on these details of the 
analysis, as long as both are carried out using the same method. We have determined m^ s and m^ s 
in pure QCD using the leading order method described here and find that the EM interactions tend 
to decrease these masses by roughly one percent (~ 0(a em )) although it is statistically less signifi- 
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cant. We are currently computing the EM splittings using PQ valence quarks and the conventional 
definition of m res to improve our determination of the quark masses. 



4. Discussion and prospects 

The systematic errors in this calculation are quickly surveyed below. In the measurements 
of neutral pions, we ignored the contribution from the disconnected diagram, which affects the 
determination of While it is possible to calculate this diagram explicitly, the clean extraction 
of a signal is likely to be difficult. 

In the study of the EM interactions, finite volume effects could be significant as the photons 
are massless. We have estimated the finite volume effects by considering the vector-saturation 
model [Q, [K], 1 1 ] and applying the physical volume of our lattices to this estimation. We find 



roughly a +10% increase in We expect a similar size of correction for the other new LECs, 
though we need to check this. 

Although it was not discussed in detail here, we also need to take into account effects of the 
third dynamical quark and non-zero lattice spacing errors, as well as effects of the quenched QED 
approximation. 
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